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PIERRE PATIE 



Abstract. The purpose of this note is to describe, in terms of a power series, the distribution function 
of the exponential functional, taken at some independent exponential time, of a spectrally negative Levy 
process § = (£t,t > 0) with unbounded variation. We also derive a Geman-Yor type formula for Asian 
options prices in a financial market driven by e^. 

Resume. Loi de la fonctionnelle exponentielle de processus de Levy assymetriques et options 
asiatiques. L'object de cette note est de donner une representation, en terme d'une serie entiere, de 
la distribution de la fonctionnelle exponentielle, consideree en un temps exponentiel independent, d'un 
processus de Levy £ spectrallement negatif, a variation infinie et pouvant etre tue. Nous en deduisons une 
formule du type Geman-Yor pour le prix des options asiatiques dans un marche financier dirige par e^. 



Version francaise abregee 

Soit £ = (£,t)t>o un processus de Levy a valeurs reelles spectralement negatif et dont les trajectoires 
sont a variation infinie. Cela signifie que £ est un processus dont les accroissements sont stationnaires et 
independants et par aillcurs lc processus n'effectue que des sauts negatifs. II est bien connu que la loi du 
processus £ est caracterisee par la loi de la variable aleatoire £ i et par consequent par l'exposant de Laplace de 
cette derniere que nous ecrivons ip. Sous les conditions H, donnees dans le corps de la note, nous proposons 
de decrire la loi de la variable aleatoire E eg = J q e^'ds, ou e q est une variable aleatoire independente de £ et 
suivant une loi exponentielle de parametre q > 0, ou nous comprenons que eo = oo. La distribution de E e(j 
apparait dans differents domaines des probabilites et egalement dans differents champs des mathematiques 
appliquees. Malheureusement, la connaissance explicite de cette loi se reduit a quelques cas particuliers, 
dont celui du mouvement brownien avec derive. A cet effet, nous indiquons l'excellent papier de Bertoin et 
Yor [5] ou une description de ces cas particuliers et des enjeux sous-jacents a l'etude de la loi de S eg sont 
dctailles. Le rcsultat principal que nous enoneons dans cette note consiste en une representation de la loi 
de S eg en terme d'une serie entiere dont les coefficients sont definis a l'aide de l'exposant de Laplace if). 
Une consequence interessante de cette representation est l'obtention d'une formule pour le prix des options 
asiatiques dans un marche financier dirige par e'. Ce resultat generalise la formule de Geman et Yor [6] 
obtenue dans le cadre du modclc de Black-Scholes. 



1. Introduction 

Let £ = (£t)t>o be a real-valued spectrally negative Levy process with unbounded variation and we denote 
its law by f y (P = Po) when £ = V & That means that £ is a process with stationary and independent 
increments having only negative jumps and its right continuous paths with left-limits are of infinite variation 
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on every compact time interval a.s. We refer to the excellent monographs [5j and [15] for background. It is 
well-known that the law of £ is characterized by its one dimensional distributions and thus by the Laplace 
exponent ip : K + — > R of the random variable £1 which admits the following Levy-Khintchine representation 

/o 
(e ur - 1 - url^ r \ <l ^)v{dr), u > 0, 
-oo 

where b G M, a > and the measure v satisfies the integrability condition J_ (1 A r 2 ) v(dr) < +oo. Since 

we excluded the case when £ has finite variation, the condition J_ (1 Ar) v(dr) < +oo is not allowed. Note, 
see e.g. [H VII, Corollary 5], that the property that £ has unbounded variation is equivalent to the following 
asymptotic 

(1.1) Km = +oo. 

u — >oo ZL 

The aim of this note is to describe the distribution function of the so-called exponential functional 

S e<i = / e is ds 
Jo 

where e g is a random variable, independent of £, which is exponentially distributed with parameter q > 0, 
where we understand e = oo. In particular, if q — 0, the strong law of large numbers for Levy processes 
gives the following equivalence 

S eo < oo a.s. -^=> E[£i] < 

and we refer to the paper of Bertoin and Yor [3] for alternative conditions, references on the topic and for 
motivations for studying the law of E e . We simply mention that this positive random variable appears in 
various fields such as diffusion processes in random environments, fragmentation and coalescence processes, 
the classical moment problems, mathematical finance and astrophysics. We also point out that the law of 
£ eij was known only for the Brownian motion with drift, see Yor's monograph [16] . and a few other isolated 
cases, see e.g. Carmona et al. [4], Gjessing and Paulsen [7] and Patie [ID]. We are now ready to summarize 
the conditions which will be in force throughout the remaining part of this note. 

H: (TTT]) holds and either q > or q = and E[£ x ] < 0. 

We mention that in [T3] the case when the condition does not hold is also considered. The remaining 
part of this Note is organized as follows. In the next Section, we state the representation of the distribution 
of E e<! in terms of a power series. In Section 3, we derive a Geman-Yor type formula for the price of Asian 
options in a spectrally negative Levy market. We end up this note by revisiting the Brownian motion with 
drift case. 



2. Main result 

Let us start by recalling some basic properties of the Laplace exponent ip, which can be found in 0. First, 
it is plain that liniu^oo ip(u) = +oo and ip is convex. Note that is always a root of the equation ip(u) = 0. 
However, in the case E[£i] < 0, this equation admits another positive root, which we denote by 8. Moreover, 
for any E[£i] € [— oo, oo), the function u i— > ip{u) is continuous and increasing on [max(#, 0), oo). Thus, it has 
a well-defined inverse function <j> : [0, oo) — > [max(0, 0), oo) which is also continuous and increasing. Finally, 
we write, for any u > and q > 0, ijj(u) — ip(u) — q, and set the following notation 

tpe(u) = ip(u + 9) and = + <K?))' 

Recalling that ip(6) = and observing that ip(4>(l)) = 0j the mappings ip$, i^M q ) are plainly Laplace exponents 
of conservative Levy processes. We also point out that i/: e (0 + ) = ip l (0 + ) > and Vw?) = ' t l ) '{ < t ) { < l)) = 
^t^j > 0. In order to present our result in a compact form, we write 

I 9 otherwise, I ipe otherwise. 
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Next, set a — 1 and a n (ip^) — (JllLi V^W) 1 , « = 1, 2, . . . In [T2], the author introduced the following 
power series 

oo 
n=0 

and showed by means of classical criteria that the mapping z i— > 2^(z) is an entire function. We refer to [12] 
for interesting analytical properties enjoyed by these power series and also for connections with well known 
special functions, such as, for instance, the modified Besscl functions, confluent hypergeometric functions and 
several generalizations of the Mittag-Lefflcr functions. To simplify the notation, we introduce the following 
definition, for any z <E C, 

O^(z) = l^z). 

Next, let G K be a Gamma random variable independent of £, with parameter n > 0. Its density is given 

— y k — 1 

by g(dy) = 6 r ^ — dy,y > 0, with T the Euler gamma function. Then, in [14] . the author suggested the 
following generalization 

0^(k;z)=E[0^ (G k z)] . 

Thus, by means of the integral representation of the Gamma function T(p) = J Q °° e~ s s p ~ 1 ds, £He(p) > 0, and 
an argument of dominated convergence, we obtain the following power scries representation 

(2.1) («; z) = —- ^(-l)"a„(^ 7 )r( K + n)z n 

' n=0 

which is easily seen to be, under the condition (jl.ip . an entire function in z. Note also, by the principle of 
analytical continuation, that the mapping p ^ 0^^{p; z) is an entire function for arbitrary zeC. We are 
now ready to state the main result of this note. 

Theorem 2.1. Assume that H holds and write S(t) — P(S e? > t),t > 0. Then, there exists a constant 
Cj > such that 

and one has, for any t > 0, 

S(t) = c^o^-t- 1 ). 

Moreover, the law of T, Bq is absolutely continuous with a density, denoted by s, given by 

s (t) = 7C 7 i- 7 - 1 a 0T (i + 7;^ 1 ), t > 0. 

We now sketch the main steps used for proving the Theorem and further details will be provided in [TJ] . 
To this end, we denote by X = {{X t )t>o, (Qx)x>o) a 1-self-similar Hunt process with values in [0, oo). It 
means that X is a right-continuous strong Markov process with quasi-left continuous trajectories and X 
enjoys the following self-similarity property: for each c > and x > 0, 

the law of the process (c _1 X ct ) t >o, under Q x , is Q x / C - 

This class of processes was introduced and studied by Lamperti [5j. In particular, Lamperti proved that 
there is a bijective correspondence between [0, oo)-valued self-similar Markov processes and real- valued Levy 
processes. Moreover, we deduce from one of the Lamperti zero-one laws that under the condition H, 
Q X (T < oo) = 1, for all x > 0, where 

T = inf{s > 0; X s = 0} 
is the absorption time of X . Then, the first key step in our proof is the following identity 

Q xt -i(H < oo) = Q x (T Q >t), x,t>Q, 

where 

(2.2) H = inf{s > 0; U s := e s X (1 _ e - s) = 0}. 
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That is H is the absorption time of U, the so-called Ornstein-Uhlenbeck process associated to X of parameter 
— 1. It is also a transient Hunt process on [0, oo). Note that the identity above is easily obtained by means 
of the self-similarity property of X and a simple time change. Observing that for t = 1, the mapping 
x i— > Q_ x (Hq < oo), defined on R + , is an increasing invariant function for the semigroup of U, we have 
thus transformed a parabolic integro-differential problem into an elliptical one. Then, specializing on the 
case when X is the self-similar Markov process associated to £ via the Lamperti bijection, we adapt, to 
the current situation, some devices which have been recently developed by the author in [T^] and [H], for 
describing the invariant function of stationary Ornstein-Uhlenbeck processes. However, several issues arise 
when dealing with the transient ones. Indeed, in the stationary case, it is an easy matter to derive some 
basic but essential properties, such as positivity and monotonicity, of the invariant functions as they are 
expressed in terms of analytical power series having only positive coefficients. As one observes from (|2.1[) . 
it is not as straightforward in the transient case and, for instance, some information about the location and 
the monotonicity of the real zeros, with respect to the argument k, of this entire function is required. This 
is achieved by combining probabilistic techniques with basic tools borrowed from complex analysis. 



3. A Geman-Yor type formula for Asian options 

In [6], Geman and Yor derived the price of the so-called Asian option in the Black-Scholes market, i.e. when 
the dynamics of the asset price is given as the exponential of a Brownian motion with drift. More specifically, 
they compute, for any K > and y € K, the following functional 

A(y,K,q)=E y [(X eq - K)+] 

in terms of the confluent hypergeometric function, in the case £ is a Brownian motion with drift. Be- 
fore stating the generalization of Geman-Yor formula, let us point out the following identity A{y 1 K, q) = 
e v A(Q, Ke~ v , q) 1 which follows readily from the translation invariance of Levy processes. We also mention 
that the fundamental theorem of asset pricing, see Delbaen and Schachermayer [5] requires that the dis- 
counted value of the asset price is a (local)-martingale under an equivalent martingale measure. However, it 
is easily checked that with the condition = r, where r > is the risk-free rate, one may carry out the 
pricing under P. We are now ready to state the following. 

Theorem 3.1. With the notation used in Theorem \2.1\ for any K > and q > tp(l), we have 
E[(S ea -if)+] = ^^ gWW ^ (rtWg )_i. J c-i ) , 

It should be emphasized that, in general, the formula above is not obtained, from Theorem 12.11 by a 
simple term by term integration. The details of the proof of this result are provided in [11] . 



4. The Brownian motion with drift case revisited 

We consider £ to be a 2-scaled Brownian motion with drift 26 £ K and killed at some independent 
exponential time of parameter q > 0, i.e. ip(u) = 2u 2 + 2bu — q and 20(g) = \j2q + b 2 — b. Note that 
V'0( g )(w) = 2m 2 + (26 + <p(q))u. Its associated self-similar process X is well known to be a Bessel process of 
index b killed at a rate q f Q X~ 2 ds. Moreover, we obtain, setting g = b + 2<f>(q), 

O- (p-x) - r (g+ 1 ) y- ( _ 1) n nP + n) 
U ^, {q) ^ X > ~ T(p) ^ { > n\T(n + g+l) [ 1 > 

= $ (p,g + 1; -x/2) 

where $ stands for the confluent hypergeometric function. We refer to Lebedev [9l Section 9] for useful 
properties of this function. Next, using the following asymptotic 

$ (p, q + 1; ~x) - r^ 1 ) -p as x _> ^ 

T{q+ 1 - p) 
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we get that — i^Trn^yu ■ Thus, we obtain, recalling that, for any q > 0, g — <j)(q) = b + <j){q) > 0, 

%(.)(*) = ^ ^f? ^" 1 ' (1 + m > e + 1; - {2trX) 
2*(9)r(^( 9 )) io ^ ; 

which is the expression [TBI (5. a) p. 105]. 

We end up this note by mentioning that in [13j some further known and new examples are detailed. In 
particular, we recover the recent result obtained by Bernyk et al. [JJ regarding the density of the maximum 
of regular spectrally positive stable processes. 
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